Time-resolved broadband analysis of slow-light propagation and superluminal 
transmission of electromagnetic waves in three-dimensional photonic crystals 
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A time-resolved analysis of the amplitude and phase of THz pulses propagating through three- 
dimensional photonic crystals is presented. Single-cycle pulses of THz radiation allow measurements 
over a wide frequency range, spanning more than an octave below, at and above the bandgap of 
strongly dispersive photonic crystals. Transmission data provide evidence for slow group velocities 
at the photonic band edges and for superluminal transmission at frequencies in the gap. Our 
experimental results are in good agreement with finite-difference-time-domain simulations. 

PACS numbers: 42.70.Qs, 78.47.+p, 41.20.Jb, 42.30.Rx 
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I. INTRODUCTION 



Electromagnetic waves propagate in straight trajecto- 
ries in homogeneous materials. This simple propagation 
changes when the wave encounters inhomogeneities and is 
scattered. Interference of the scattered waves can modify 
significantly their propagation. The careful engineering 
of the inhomogeneities may lead to the inhibition of the 
propagation in certain directions and to the enhancement 
in others. It is also possible to speed-up P, S S 0> IE S 
or slow-down 0, BM ITll lTp the propagation of the 
wave and even to localize it [l^l- The full control of the 
propagation of electromagnetic waves is intensively pur- 
sued since this control will lead to new conce pts in nu- 
merous fields such as information processing hj, laser 
physics I bio-sensing |0| and quantum optics Tt'I . 

One of the most prominent candidates for the control 
of the propagation of electromagnetic waves are photonic 
crystals. Photonic crystals are periodic structures of two 
or more materials with different refractive indices [islIT^ . 
Interference of waves scattered at different lattice planes 
of the crystal determine its optical properties. Depend- 
ing on the structure and if the scattering is strong enough 
a photonic band gap (PBG) might be created. A PBG 
is a frequency range in which no optical modes are al- 
lowed and, consequently, the propagation of waves in 
this range is forbidden. Particulary attractive is the 
extraordinary optical dispersion that photonic crystals 
exhibit, which leads to a pronounced variation of the 
wave's group velocity. The resonant scattering at the 
band edges of the gap gives rises to a strong reduction 
of the group velocity and a significant group velocity dis- 
persion 0. This slow propagation, which increases the 
interaction time between the wave and the crystal, con- 
stitutes the basis for the development of sensors for gases 
and bio-molecules and more sensitive non-linear com- 
ponents ■ A wave with a frequency in the PBG inci- 
dent on a photonic crystal is reflected. However, if the 
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crystal has a finite thickness, a fraction of the wave's 
intensity is transmitted. While the transmitted inten- 
sity decreases exponentially with the crystal thickness, 
the transmission time is independent of this thickness, 
which may lead to superluminal transmission. Tunneling 
processes through photonic structures are ofOTeat scien- 
tific and technological interest 0, [2^ I3 . Among 
the reasons driving this interest are the similarities with 
the tunneling of quantum particles and the ongoing in- 
tense discussions on the propagation in the superluminal 
regime 23, 24]. 

In spite of the fundamental and technological relevance 
of dynamic properties of electromagnetic waves in pho- 
tonic crystals most of the investigations so far focus only 
on stationary measurements. This focus is mainly im- 
posed by the complexity of phase sensitive techniques at 
high frequencies. Reports on the group velocity of waves 
in optical photonic crystals are thus seldom and limited 
to narrow frequency rang es 0] , to weakly scattering pho- 
tonic crystals SHISllfl, or to low dimensional photonic 
structures HHIIMBQ 

In this article we present a broadband time domain in- 
vestigation of the propagation of electromagnetic waves 
in strongly-scattering-three-dimensional photonic crys- 
tals of variable size. In our investigation we use terahertz 
time-domain spectroscopy, a technique capable to gener- 
ate ultrashort broadband pulses of THz radiation and 
to detect the amplitude and the phase of these pulses. 
Through the direct access of amplitude and phase an un- 
precedented precision in the analysis of the propagation 
of electromagnetic waves in photonic crystals is achieved. 
The extinction coefficient is obtained from the exponen- 
tial decrease of the transmitted amplitude with the thick- 
ness of the crystal. From the phase analysis of the trans- 
mitted wave the wave number, the effective refractive 
index and the group velocity are derived. Very strong 
dispersion at the edges of the gap and anomalous disper- 
sion at frequencies in the gap are observed. Even in a 
photonic crystal with a thickness of only four unit cells 
the group velocity is dramatically reduced at the band 
edges by almost a factor 15 with respect to the speed of 
light in vacuum, while anomalous dispersion in the gap 
photigifttflgaagts^ to superluminal transmission. 
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II. SAMPLES 

The photonic crystals investigated are formed by in- 
trinsic sihcon and air. Intrinsic siUcon has a very large 
refractive index n ~ 3.4, and virtually no absorption at 
terahertz frequencies. Photonic crystals are fabricated 
by micromachining thin wafers of silicon, and piling them 
to form a layer-by-layer structure [23, Hg] • Such a struc- 
ture consists of a transverse stacking of arrays of dielec- 
tric rods with a face-centered-tetragonal lattice symme- 
try [27! l28j . This type of cr ysta l structure exhibits a 
very large photonic band gap [23 ■ The micromachining 
of our crystal was done with a programmable dicing saw, 
as described in more detail in Ref. 26]: parallel grooves 
spaced by 185 /xm were diced on one side of a wafer with 
a thickness of 138 ± 3 /j,m. The thickness of the grooves 
is 110 fim and the depth one half of the thickness of the 
wafer, i.e., 69 ^m. Perpendicular grooves to these with 
the same thickness and depth were diced in the opposite 
side. This method produces a square grid of apertures. 
To form the photonic crystal structure several of these 
grids are piled. Consecutive grids are shifted by one half 
of the period and two of them form a single unit cell in 
the stacking direction. The dimension of this unit cell 
is thus 276 /im. The calculation of the band structure 
of this specific photonic crystal predicts a PBG centered 
at 0.5 THz. Figures n^a) and (b) are scanning electron 
microscope images of a layer-by-layer structure with a 
thickness of 4 unit cells. Figurema) is a side view of the 
structure along the stacking direction, while (b) is a top 
view of the same photonic crystal. 



III. TERAHERTZ TRANSMISSION THROUGH 
PHOTONIC CRYSTALS 

A. experimental setup 

For the experiments we use a THz time-domain spec- 
trometer [291 . This setup is schematically depicted 
in Fig. 13 A train of femtosecond pulses from a 
Ti:sapphire laser is split to generate and detect single- 
cycle terahertz pulses. One beam, the so-called pump 
beam, is incident on an InGaAs surface field emitter. 
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FIG. 2: Schematic representation of the experimental setup. 
Terahertz pulses are generated on an InGaAs surface with a 
laser pulse. These pulses are coUimated and incident normal 
to the crystal's surface. The transmitted amplitude is focused 
onto a photoconductive antenna which is gated by a second 
laser pulse. By varying the path length difference between the 
two optical paths the THz amplitude is detected as a function 
of time. To avoid water vapor absorption the THz setup is 
enclosed in a box which is continuously fluxed with nitrogen. 

generating THz radiation. The THz pulses are collected 
and coUimated by a parabolic mirror placed at the fo- 
cal distance to the emitter. This THz beam is incident 
normal to the surface of the crystal, and the transmit- 
ted amplitude is focused with a parabolic mirror onto a 
photoconductive antenna. The antenna is gated with the 
pulses from the second beam of the Ti:sapphire laser. By 
means of a computer controlled delay stage, the length 
difference between the generation and detection optical 
paths is varied and the THz field amplitude is detected as 
a function of time with sub-picosecond time resolution. 
Water vapor absorbs THz radiation, therefore, the THz 
setup is enclosed in a box which is continuously fluxed 
with nitrogen gas. 

Even though the measurements were done at only one 
direction, it has been demonstrated by Imhof et al. that 
ID models are not able to explain the measurements of 
the group velocity near the band edges of 3D photonic 
crystals 



(a) (b) 




FIG. 1: Scanning electron microscope photographs of a Si 
layer-by-layer photonic crystal with a thickness of 4 unit cells. 
Photograph (a) is a side view, while (b) is an upper view. 



B. Transmission measurements 

We have measured the transmission of THz pulses 
through photonic crystals with different thickness rang- 
ing from 1 to 4 unit cells. The experimental transients of 
the transmitted field amplitude are plotted in Fig. |2Ia). 
The reference pulse, i.e., the setup response, is also plot- 
ted in the same figure. 

The measurements are compared with numerical sim- 
ulations using the finite-difference-time-domain FDTD 
method. These simulations are plotted in Fig.|2Ib), and 
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FIG. 3: Terahertz pulses transmitted through photonic crys- 
tals with different thickness. The thickness is indicated in the 
figure as the number of unit cells (1 unit cell = 27Q^,m). Mea- 
surements are plotted in (a), while the corresponding FDTD 
simulations are displayed in (b). For clarity a vertical offset 
has been introduced to the different pulses. The uppermost 
pulses are the reference, i.e., the instrumental response of our 
setup. 



were performed with a commercial software (FullWave by 
RSoft). The spatial resolution of the simulations is 3 /im 
and the temporal resolution is 5 fs. The incident field was 
approximated to a plane wave by using a Gaussian field 
distribution with a width much larger than the size of 
the simulation region. The temporal shape of the input 
THz field was a single-cycle pulse, with a center wave- 
length of 500 /im. The transmitted field was calculated 
at a distance of 2.3 mm from the source. 

As can be seen in Fig. O the transmitted pulse experi- 
ences a significant dispersion growing with the thickness 
of the crystal. Signatures of the PBG formation become 
directly visible in the time domain data. Particulary in 
the thickest crystals, low frequencies, belonging to the 
lower photonic crystal band, and high frequencies, asso- 
ciated to higher bands, are separated in time domain. 

To analyze the transmission we calculate the com- 
plex Fourier transform of the time-domain transmission 
through the crystal and normalized it by the reference 
measured without the sample. This transformation leads 
to a transmission T{v) of the form 



(1) 



The modulus t{v) is the spectrally resolved field trans- 
mission coefficient, and the argument A0(i/) is the phase 
difference between the transmission through the crystal 
and the reference. Equation can be written as 



where fcos is an effective propagation constant for the 
transmission through the crystal and L is the thickness 
of the crystal. The effective propagation constant is given 
by 



k,s{v) = + ikM 
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where ricff is the effective refractive index and k^s the 
extinction coefficient of the crystal. The transmission is 
thus proportional to 



(4) 



The amplitude transmitted through the crystal decreases 
exponentially with the crystal thickness 



while the phase difference is given by 
27r 



A 



(noff(i^) -l)L 



(5) 



(6) 



ft is important to mention that the transmission is mea- 
sured outside the photonic crystal. Therefore the prop- 
agation constant has a real component or wave number 
k'^ff{v) even at frequencies within the gap. In section lTlI (]l 
the analysis of the phase and the dispersion of the pho- 
tonic crystals are discussed in detail. 

The amplitude transmission coefhcient t and the phase 
difference of a photonic crystal with a thickness of 3 
unit cells are plotted in Fig. 0] The circles in this figure 
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FIG. 4: (a) phase difference between the transmission through 
a photonic crystal with a thickness of 3 unit cells and a ref- 
erence measurement, (b) Amplitude transmission spectrum 
normalized by the reference. The wave incidences normally 
to the crystal surface. The solid lines correspond to a FDTD 
simulation. 
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FIG. 5: Amplitude transmission at the center frequency of 
the gap as a function of the thickness of the crystaL The 
sohd hue is an exponential fit to the measurements. 



correspond to the measurements and the solid lines are 
results of the FDTD simulation. All the crystal param- 
eters in the simulation, i.e., structure, dimensions and 
refractive index, are known, and the remarkably good 
agreement between simulation and measurement is ob- 
tained without adjusting any parameter. The transmis- 
sion (Fig. El^b)) drops strongly in the frequency range 
defined by the gap, i.e., 0.4-0.6 THz. At these frequen- 
cies the phase difference (Fig. EJa)) is nearly constant. 
The interference pattern that is apparent by the oscilla- 
tions in the transmission at low and high frequencies are 
Fabry-Pcrot resonances due to multiple reflections at the 
crystal interfaces. 

The extinction coefficient KcS describes the attenua- 
tion of the wave amplitude as it propagates through the 
medium. In photonic crystals made of non-absorbing 
constituents, such as those investigated here, the extinc- 
tion of the transmitted amplitude at frequencies within 
the gap is due to destructive interference of multiply scat- 
tered waves. At these frequencies the wave is Bragg re- 
flected and only a small fraction is transmitted through 
the crystal. According to Eq. ((SJ the amplitude decreases 
exponentially with the crystal thickness and with an at- 
tenuation length La = 2tt^ ff • This attenuation length 
is a direct measure of the pliotonic strength of the crys- 
tal. Weakly scattering photonic crystals have attenua- 
tion lengths of the order of tens of unit cells, while La in 
strongly scattering crystals is on the order of a few unit 
cells. 

Figure ]E\ depicts the minimum transmission centered 
at the gap frequency of 0.5 THz as a function of the 
thickness of the crystal. The line in Fig. |31 is an expo- 
nential fit to the measurements from which we obtain an 
amplitude attenuation length of 195 ±20 fj,m. This at- 
tenuation length is even smaller than one unit cell (~ 0.7 
times the lattice constant), indicating an extraordinarily 
strong scattering in Si layer-by-layer structures fs^, Wl] ■ 

The extinction coefficient is represented in Fig. |^ for 
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FIG. 6: effective extinction coefficient of Si photonic crystals. 



different frequencies. Radiation with a frequency outside 
the gap can propagate through the crystal being only 
weakly reflected at the interfaces. For these frequencies 
we have k^s — 0. At frequencies within the gap the trans- 
mitted amplitude decreases exponentially with L leading 
to a positive KcS- The largest value of Kcff, or equiva- 
lently the smallest La, is at the central frequency of the 
gap. 



C. Phase analysis 

The dispersion properties of the electromagnetic wave 
propagating through the crystal were obtained as de- 
scribed in Refs. |32l and in the previous section. 
The measured phase difference between the transmission 
through the photonic crystal and the reference was used 
to calculate the wave number 
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where L is the thickness of the photonic crystal, Co is 
the speed of light in vacuum and uj = 2'itv the angular 
frequency. The measured wave number of a layer-by-layer 
structure of 3 unit cells is shown in Fig. [3 with circles. 
The solid line in this figure represents the data from the 
FDTD simulations. The difference between measurement 
and simulation at high frequencies can be ascribed to 
small imperfections and misalignments in the structure. 

The effective refractive index of the crystal ricff, de- 
fined as ncff(z^) — i^A;'(t/), is plotted in Fig. |H|versus the 
frequency. The FDTD simulation is represented in the 
same figure with a solid line. There is a distinct increase 
of rics at the low frequency band edge, i.e., 0.4 THz, fol- 
lowed by a remarkable decrease in the gap. At the high 
frequency band edge, i.e., 0.6 THz, n^s rises again. As 
we are going to see, the increase of the effective refractive 
index at the band edges gives rise to a strong reduction 
of the group velocity, while the anomalous dispersion in 
the gap leads to superluminal transmission. 
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FIG. 7: Frequency versus wave number along the stacking 
direction of a Si photonic crystal with a thickness of 3 unit 
cells. The solid line is a FDTD simulation. 
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FIG. 9: Inverse of the group velocity multiplied by the speed 
of light in vacuum versus the frequency of a Si photonic crystal 
with a thickness of 3 unit cells. The dashed line is a guide to 
the eye. The solid line is a FDTD simulation. 



To calculate the group velocity Vg, we use the standard 



definition 



where 



dAct> 



is the group time 



delay. With the expressions of k' and Uoff, the effective 
group velocity can be written in the familiar form 



Co 



rieff + LodueS /duj 



(8) 



The dispersion in Vg introduced by the photonic crystal 
is described by the second term of the denominator, i.e., 
the derivative of rics with respect to the frequency. 

The inverse of the group velocity multiplied by Co, i.e., 
the denominator of Eq.|Sl is plotted in Fig. with circles, 
and the corresponding FDTD simulation with a solid line. 
The large dispersion at the band edges due to the reso- 
nant scattering gives rise to a strong reduction of the 
group velocity. The oscillations of the group velocity at 
low and high frequencies that are clearly visible in the 
simulation are due to the multiple reflections at the crys- 
tal interfaces. The anomalous dispersion at frequencies 
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FIG. 8: Frequency dependence of the effective refractive index 
of a Si photonic crystal with a thickness of 3 unit cells. The 
solid line is a FDTD simulation. 



in the gap, which is characterized by the negative value 
of ducs/duj, increases the group velocity as defined by 
Eq. IHl If the crystal is thick enough this group velocity 
may be larger than the speed of light in vacuum. This 
result can be appreciated in Fig. IHl where the average 
value of Co / Vg in the gap of a crystal with a thickness of 
only 3 unit cells is 0.4 ± 0.8. 

Although the analysis presented above is the most 
accurate for the study of the dispersion characteris- 
tics of photonic structures, time- frequency representation 
methods are more intuitive and might be useful in this 
study ^3^1 . As an example of these methods, we plot in 
Fig.^lthe short time Fourier transform STFT ,32j of the 
transmission through the photonic crystal with a thick- 
ness of 2 unit cells. The STFT is calculated by applying 
a window or a band pass filter to the time domain signal 
and Fourier transforming the signal within this window 
to obtain the spectral contain for a certain time delay. 
The process is repeated by shifting the time window over 
the entire time interval of the measurement. In this way 
temporal and spectral information are simultaneously re- 
trieved. The contour plot of Fig. llUf al is the transmit- 
ted amplitude through the photonic crystal, where the 
blue stands for no transmission and red for high trans- 
mission, as a function of time and frequency, calculated 
with a Hanning time window 8 ps wide. The THz tran- 
sients are plotted in Fig. Iior b). where the red curve is 
the transmission through the crystal and the blue line 
corresponds to the reference. The amplitude transmis- 
sion coefficient is plotted in Fig. llOf c) . Consistently with 
the phase analysis, the time delay is the largest at the 
band edges, where the group velocity reaches its mini- 
mum value. At frequencies within the gap the transmis- 
sion decreases significantly and the time delay becomes 
smaller. 

Figure [TTI summarizes the salient features of the group 
velocity and presents its dependence on the thickness of 
the crystal. Fig. ^] (a) displays the inverse of Vg mul- 
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FIG. 10: (a) Short time Fourier transform of the THz trans- 
mission through a Si photonic crystal with a thickness of 2 
unit cells. The transmitted pulse and the reference are repre- 
sented in figure (b) with a blue and red lines respectively. The 
amplitude transmitted through the crystal normalized by the 
reference is plotted in figure (c) as a function of the frequency. 



tiplied by Cq at the high frequency band edge, and (b) 
represents the average value of Co/fg in the gap. The 
circles in these figures correspond to the measurements, 
while the open triangles to the FDTD simulations. The 
group velocity at the band edge is rapidly reduced as the 
thickness of the crystal increases. With only four unit 
cells this group velocity is almost a factor 15 slower than 
the speed of light in vacuum. On the other hand, the 
group velocity in the gap rapidly increases as the thick- 
ness of the crystal is augmented. For a photonic crystal 
of 3 unit cells the group velocity is roughly two times 
larger than the speed of light in vacuum. The transmis- 
sion at gap frequencies of the thickest photonic crystal 
was lower than the noise level of our setup. For an infi- 
nite and defect-free photonic crystal one should expect a 
zero group velocity at the band edges and a no transmis- 
sion in the gap. 



through strongly scattering photonic crystals with dif- 
ferent thickness over a large bandwidth covering more 
than an octave below, at and above the bandgap. At fre- 
quencies within the gap we observe an attenuation length 
smaller than o ne unit cell and anomalous dispersion. The 
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FIG. 11: Inverse of the group velocity multiplied by the speed 
of light in vacuum as a function of the thickness of the pho- 
tonic crystal in unit cells. The circles are measurements, while 
the open triangles correspond to FDTD simulations, (a) cor- 
responds to the high frequency edge of the gap, while the data 
in (b) are at frequencies in the gap. The dashed line in this 
figure indicates the speed of light in vacuum. The solid lines 
are guides to the eye. 



strong dispersion introduced by the crystal gives rise to 
extremely low group velocities at the band edges and to 
superluminal transmission in the gap. 
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